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We present an evaluation of the tt cross section near threshold at next-to-next-to-leading logarith- 
mic accuracy, using a two-step matching procedure. QED corrections are taken into account as 
well and are shown to be numerically important. Finally we give an outlook on how to further 
improve the theoretical predictions with particular emphasis on how to consistently include finite 
width effects. 
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1. Introduction 

A future International Linear Collider (ILC) offers the opportunity to study the top quark with 
unprecedented accuracy. To fully exploit the potential of an ILC in this respect, it is essential that 
a dedicated measurement of the cross section for the production of a top antitop quark pair close 
to threshold is made. Such a threshold scan allows for an extremely precise measurement of the 
top-quark mass and yields information on the top-quark width and the top-Higgs Yukawa coupling. 

In order to take full advantage of future experimental data, theoretical computations have to 
match the predicted experimental precision. As is well known [|l|], a perturbative expansion in 
the strong coupling a^ alone is not adequate due to the presence of terms {ajv)" at ?i-th order in 
perturbation theory, where v <C 1 is the (small) velocity of the top quarks near threshold. These 
terms have to be resummed and the theoretical computations are organized as a double expansion 
in the two small parameters a^ ~ v. Thus, a next-to-leading order (NLO) result includes all terms 
that are suppressed by one power of either «« or v, whereas a next-to-next-to-leading order (NNLO) 
result contains all terms suppressed by either aj, a^ v or v^. 

Such an expansion is done most efficiently by using an effective theory approach. From a tech- 
nical point of view this amounts to using the threshold expansion [^ to split the Feynman diagrams 
into contributions due to various modes and then integrating out unwanted modes. The modes to 
be considered are hard modes (with momentum of the order of the heavy quark mass k^ ~ m,), soft 
modes (with momentum of the order of the typical momentum of the non-relativistic top quarks 
k^ ^nitv), ultrasoft modes (with momentum of the order of the typical energy of the non-relativistic 
top quarks k^ ^ nit v^) and, finally, potential modes (with k^ ^ nifV^ and k r^ rritv). In a first step hard 
modes are integrated out, resulting in non-relativistic QCD (NRQCD) ^. NRQCD is matched to 
QCD at a hard scale /X/, ^ mt. In a second step, soft modes and potential gluons (and massless 
quarks) are integrated out. The resulting theory is called potential NRQCD (pNRQCD) [Q] and is 
matched to NRQCD at a soft scale pis ~ rut v. This theory consists of heavy quarks with energy 
E ~ m,v^ and momentum k ~ m,v, interacting through potentials (to be interpreted as matching 
coefficients) and still dynamical ultrasoft gluons with momentum k^ ~ /x,„ r^ nitV^. 

This approach has been used by several groups to compute tt-paii production near threshold 
at NNLO (for a review see [^]). If the cross section is expressed in terms of a threshold mass [^, ^ 
rather than the pole mass, the position of its peak is very stable and can be predicted with a small 
theoretical error. This will allow to determine the top threshold mass and ultimately the top MS- 
mass with a very small error. The situation is much less favourable regarding the normalization 
of the cross section. The corrections are huge and the scale dependence at NNLO is larger than at 
NLO, indicating that this quantity is not well under control at NNLO. 

The bad behaviour of the fixed order perturbation theory is due to the presence of large log- 
arithms log V. In order to improve the situation, these logarithms have to be resummed. Counting 
a log V ^ 1 in a NLO and NNLO calculation produces a result of next-to-leading logarithmic (NLL) 
and next-next-to-leading logarithmic (NNLL) accuracy respectively. At leading order, there are no 
logarithms. Thus the leading order (LO) result coincides with the result at leading logarithmic ac- 
curacy (LL). Within the effective theory approach the logarithms are resummed by computing the 
anomalous dimensions of the matching coefficients and solving the corresponding renormalization 
group equations. 
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The renormalization-group improved (RGI) cross section for ff pair production close to thresh- 
old has been computed [§] a few years ago using vNRQCD [^, an approach that is somewhat 
different from the two step matching procedure QCD —>■ NRQCD —>■ pNRQCD mentioned above. 
In vNRQCD the relation between the soft and ultrasoft scale is fixed from the beginning by setting 
IJ.US = IJ'j/nit and there is only one step in the matching procedure, QCD -^ vNRQCD. These com- 
putations showed that the log v terms are numerically very important and improve the behaviour of 
the perturbative series considerably. 

So far no calculation of tt pair production at NNLL accuracy using a pNRQCD approach has 
been done, even though most partial results were available for quite a while. In Section ^ we 
present such a computation []T0[]. In addition we include QED corrections which turn out to be 
numerically relevant. The most important missing terms are identified and prospects for future 
improvements are given in Section ^ A particularly important problem is the inclusion of effects 
due to the instability of the top quarks. This will be discussed in Section ^ before the conclusions 
are presented. 

2. Current status: tt at NNLL 

The starting point for the RGI evaluation of the threshold cross section is the fixed-order 
calculation presented in [pi]]. The cross section obtains contributions from 7 and Z exchange. 
In order to simplify the discussion we ignore the Z exchange in what follows. The cross section at 
center of mass energy ^/s can then be written as 

a^s) = ^R{s) (2.1) 

where a is the electromagnetic coupling and the R ratio is expressed as a correlator of two heavy- 
quark vector currents j'^ix) = Qy^Q{x) 



2 



R{s) = ^Im 



d\e'^'^{0\T{f{x)i^{0)m 



(2.2) 



with q^ = s = {2mt +eY and Ct the electric charge of the top. The current is expressed in terms of 
two-component spinor fields i/a^ and % resulting in j ~ ci;^^a'VA — C2/(6m,^);^^a' (/D) i/a+ .... 
The contributions of the hard modes, which are to be integrated out first, go into the coefficient 
functions di{ji.h) of the NRQCD Lagrangian and into the matching coefficients of the current 
Cid^h)- For most coefficients a leading-order result is sufficient, however, ci has to be computed to 
NNLO [|l2|]. Counting a ~ aj, the one-loop exchange of a hard photon contributes at NNLO and 
is taken into account by 

ci(M/,)-ci(At,)-^^^^^ (2.3) 

Once the hard modes are integrated out, we are left with soft, potential and ultrasoft modes. How- 
ever, near threshold the only allowed external states are potential heavy quarks and ultrasoft gluons. 
Integrating out the remaining modes results in the pNRCQD Lagrangian, which has the structure 

^pNRQCD = ^pSrqcd + Jd'r [w^r' y] {x + f) 5V [X^T'^X] (x) (2.4) 
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where the leading order Lagrangian 

^pSIqcd = '^'(^■^" + ^ + ^T')'/^ + ^'('^''-^ + ^T')^ (2-5) 

+ Jd'r [y'T^y] (x + f) (-y) [X'^T"X] (x) 

includes the leading Coulomb interaction and takes into account the effects due to the finite width 
of the top quark by the shift £ ^ £" + /F,. It is known that this does correctly take into account the 
width effects at leading order, but not at NNLO, an issue we will come back to in Section ^. 

The starting point is the free Green function of the Lagrangian -^pnrqcd which, from a dia- 
gramatical point of view, includes all ladder diagrams with the exchange of an arbitrary number of 
potential gluons. The imaginary part of the Green function at the origin is related to the R ratio by 

R^E) = ^1^ (^c?-ciC23^) ImG(0,0,£) (2.6) 

where A'c = 3 is the color factor. In order to regularize the ultraviolet singularity present in the 
Green function, we perform the calculation using dimensional regularization inD = d+l =A — 2e 



dimensions [11]. After minimal subtraction we obtain [13] 



''-("•"•'^'"/iSF^'^'*^""'"'^'''' 



asCfin} /I 1 -4mtE 1 . ^^ i 



471 V2A 2 n^ 2 



where A = Cf «« / {2y/—E/mt ) with Cf = 4/3. Once the potential, 5V, is known to the desired 
accuracy, the NNLO corrections are calculated as single or double insertions 

J (2^(2^(2^(2^^^'(^^'^^'^)^^(^^-'^^)^'-(^i'^^'^) ^'-'^ 

Such insertions can result in divergences. Since the whole calculation is performed using dimen- 
sional regularization, it is important that the potentials, to be interpreted as matching coefficients, 
are computed consistently in D dimensions. Using the result for the two-loop static potential [|l4|]. 



the potential has been computed in Ref. [|rT|]. Here we include an additional contribution in 5V, 
the electromagnetic Coulomb term 

5V^5V-'-^ (2.9) 

This will give rise to a NLO term ^ a/v from single potential photon exchange and a NNLO 
term from double potential photon exchange. We also remark that in Eq. ( p!4| ) the leading ultrasoft 
interactions are included [ [131 ] even though they only contribute beyond NNLO. They are given by 
the chromoelectric dipole operator x-E{t,0) where the electric field is understood to be multipole 
expanded. An additional term containing A°(?,0) is not displayed since it can be gauged away. 

We now turn to the discussion on how to resum the logarithms in this approach. After matching 
NRQCD to QCD at the hard scale jLt/,, the renormalization group equations are used to evolve 
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the matching coefficients di from di{}JLh) to the soft matching scale di{}JLs), thereby resumming 
logarithms of the form logpLh/pLs- The matching coefficients of the single heavy quark sector can 
be taken from heavy quark effective theory [jTsI]. The RGI coefficients of the four-heavy-quark 
operators have been computed in Ref. [jl^]. 

Matching pNRQCD to NRCQD at the soft scale }Xs then results in potentials that depend on 
as{lJ.s) and the matching coefficients di{}i.s)- Again, renormalization group equations are used to 
evolve the potentials down to ultrasoft scales pL^s, resumming log Us / l-ius [[161] . This procedure is 
complicated by the fact that higher dimension operators of the single heavy quark sector can mix 
through potential loops into lower dimension operators of the heavy quark-antiquark sector. Thus, 
to obtain the NLL (NNLL) running of ci the NLL (NNLL) running of (some of the) single heavy 
quark operators has to be taken into account. While all other coefficients are known to an accuracy 



sufficient for a NNLL calculation, ci is only known to NLL []170. Some NNLL terms of ci have 



been computed 018| ] but we stress that strictly speaking the term NNLL that is commonly used to 
describe the accuracy of the result is not valid. 

From the discussion in the introduction it is clear that the hard, soft and ultrasoft scales are 
correlated /x,u ~ l^j/l^-h- As mentioned above, in the vNRQCD approach the scales are fixed by 
Hus = pL^ /nit. While the two-step matching procedure in principle does allow to relax this condition, 
in the current approach, the relation }Xus = l^j/lJ-h has been used in solving the renormahzation 
group equations. Thus, it is not possible to vary Hus, i^k and /x^ independently around their natural 
values. This is somewhat unfortunate, as it would be preferable to have a formulation where the 
correlation of the scales is not rigid and which would allow to treat and vary the various matching 
scales independently. This would allow to obtain a more realistic indication of the theoretical error, 
in particular since some ultrasoft logarithms at NNLL are missing in the present results. For the 
moment we just have to keep this in mind when trying to assign a theoretical error. 




H,= AO 




1 2 

-4^-2 mps 



Figure 1: Threshold scan for tt using nit^ps = 175 GeV. The left figure shows the fixed order results, LO, 
NLO and NNLO, whereas the in right figure the RGI results LL, NLL and NNLL are displayed. The soft 
scale is vai-ied from /X.(=40 GeV to ^i=80 GeV 

In Figure |l| the fixed-order results are compared to the RGI results |[l^] obtained using the 
procedure described above and using the PS mass [^. In order to obtain a first rough estimate 
of the theoretical uncertainty we show the cross sections as bands obtained by variation of the 
soft scale in the region 40 GeV < /i^ < 80 GeV and setting /X/, = m,. As was to be expected 
from previous results [^], the scale dependence is much reduced once the logarithms are taken into 
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account. We also note that the size of the corrections decreases for the RGI results, in particular 
the NNLL band is much closer to the NNL band, and the scale dependence reduces from LL to 
NLL to NNLL. However, the NNLL band does not overlap with the NLL band, indicating that a 
theoretical error estimate relying on the scale variation in this plot alone is too optimistic. 

In the fixed-order calculations the by far dominant uncertainty came from the variation of the 
soft scale. Since this problem is much less severe after the resummation of the logarithms we 
have to be more careful with other sources of uncertainties. In particular, the missing ultrasoft 
contributions make it important to consider the dependence on the other scales as well. Since 
/i„.v is related to jj.^ by Hus = l^j/l^h, we consider in the left panel of Figure ^ the dependence 
on /I/, for Us = 40 GeV. Variation of the hard scale around its natural value /i/, = nit by choosing 
100 GeV </!/,< 250 GeV results in a scale dependence that is considerably larger than the soft 
scale dependence. It is to be expected that the situation improves once all ultrasoft logarithms at 
NNLL are taken into account, but at this stage the rather large dependence of the cross section on 
Hh has to be taken into account if a theoretical error is assigned. We also note that the dependence 
on fXh only enters at NLL, thus the LL "band" in the left panel of Figure]^ is simply a line. 



^lh= 100 -250 GeV 



0.4 





Eps=Vr-2r 



Eps=v's"-2r 



Figure 2: Left panel: Variation of the hard scale Hf, in the threshold scan for m,ps — 175 GeV. The hard 
scale is varied from 100 GeV < jj.;, < 250 GeV. Right panel: the effect of including QED corrections at NLL 
and NNLL. There are no QED corrections at LL. 

Finally we turn to the right panel of Figure ^ where we display the effects due to the QED 
corrections. As mentioned, they enter at NLL (thus there is no effect on the LL curve) and they are 
large. They change the normalization by up to 10% and result in a shift in the extracted MS mass 
of up to 100 MeV, making their inclusion imperative. 

3. Beyond NNLL 



The theoretical error on the peak of the cross section and, with the current RGI results, also its 
normalization is in rather good shape. Depending on the details of the error analysis, the theoretical 
error in the determination of the top mass is in the range of 50-150 MeV. The uncertainty in the 
normalization of the cross section has been quoted as 6% in Ref [^ whereas Figure ^ seems to 
indicate a somewhat larger value. In any case it is desirable to further reduce the theoretical error, 
in particular in the determination of the normalization of the cross section, in order to obtain the 
best possible measurements of the top width and maybe even the top Yukawa coupling. 
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An obvious first step is to complete the NNLL determination of the current matching coeffi- 
cient ci- Once this is completed, there are two main directions for further improvements. The first 
concerns the correct implementation of finite width effects and will be discussed in Section ^ The 
other is the inclusion of connections of even higher order. 

A full NNNLO calculation requires the evaluation of the matching coefficients to the desired 
order, as well as the computation of the insertions of the potentials. Several partial results have 
already been obtained. For a start, the logarithmically enhanced NNNLO terms have been obtained 
a while ago [p^. These terms are included in the NNLL result. Furthermore, the matching coeffi- 
cient for the non-analytic 1/r^ potential has been computed to the required accuracy [20] and the 
NNNLO contribution due to multiple insertions of Coulomb potentials has been computed [p|]. 
However, the calculation of the static potential and the matching coefficient ci of the current are 
major obstacles on the way to a full NNNLO evaluation. 

Of particular importance are the ultrasoft effects. As mentioned before, they enter only at 
NNNLO and arise due to the emission and absorption of an ultrasoft gluon through the chromoelec- 
tric dipole operator. Thus one of the couplings is a{jJLus)- Since the ultrasoft scale is rather small, 
it can be argued that these contributions are potentially particularly important NNNLO terms and 
their evaluation would be an important partial result. Note that this contribution has an ultraviolet 
singularity, resulting in logarithmically enhanced NNLL terms [p^]. 

The importance of NNNLO contributions can also be seen from the fact that multiple insertions 
of the Coulomb potential are numerically significant. In Figure [T] the soft scale has not been allowed 
to run below 40 GeV. The reason is that for /ij < 30 GeV the scale dependence get much worse. 
This bad behaviour is due to the Coulomb potential and it has been shown that including multiple 
insertions of the Coulomb potential the situation can be rectified [^l]]. 

4. Finite width effects 

The width of the top quark, Y, is a leading-order effect that has to be taken into account in the 
propagator of the non-relativistic heavy quarks. This can be seen by noting that the non-relativistic 
propagator [E — p^/(2m,)]^' ~ [m^v^]^' scales in the same way as the width Yt ~ nitaew ~ /m,v^. 
Thus the suppression due to the insertion of the self energy is compensated by the enhancement of 
an additional non-relativistic propagator and the corresponding resummation yields 

1 _/ -iY, Y 1 

^' ' ^ (4.1) 



E-py{2mt)^\E-py{2m,)J E + iY, - p^ / {2m,) 

This is the justification at leading order of the replacement E ^ E + iYf The question is how to 
go about to include the finite width effects beyond leading order. This question is actually relevant 
for a large number of processes involving unstable particles, but it is particularly pressing here, 
given the unprecedented precision required in the theoretical predictions. We recall that the error 
of the top quark mass 8mt ~ 100 MeV, thus 5m, <^ Y,. Therefore, a systematic approach on how 
to include finite width effects is mandatory. 

To start with we have to emphasise that the very notion of a tt cross section breaks down. 
Taking into account the width effects means we have to take into account the electroweak decay 
of the top quark. Thus we are dealing with the process y* -^ W^^b W^b. Even if we neglect 
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the decay of the W bosons, this opens up the possibility of many additional radiative correction 
processes. There are gluons connecting the intermediate top quarks with the decay products (non- 
factorizable corrections), there are genuine electroweak corrections and, on top of this all, there are 
radiative corrections due to photon exchange from any of the decay products with any other charged 
particle, even the incoming electrons. Some of the electroweak corrections have been taken into 
account through absorptive parts in the NRQCD matching conditions [^] and have been shown 
to be numerically relevant. However this does not take all effects into account. In particular QED 
radiative corrections which are an order a effect and thus contribute at NNLO have to be considered 
as well in a complete solution to the problem. 

The key to the solution is the observation that the situation regarding the finite width correc- 
tions is actually very similar to NRQCD. In QCD we have at each order n of perturbation theory 
terms of the form (as/v)". Thus the suppression by the coupling is compensated by a kinematic 
enhancement. This is exactly the same as in Eq. (p^. The Coulomb singularity v = has its 
analogy in the propagator pole of an on-shell E = p'^/{2mt) stable Ff = top quark. Resummation 
of the {as/v)" terms results in a function that is well defined at v = in the same way as the right 
hand side of Eq. ( p~I| ) is well defined for E = p^/{2mt). In NRQCD, the problem of the Coulomb 
singularity is solved by systematically expanding in both small parameters a^ and v. Thus the basic 



idea is to do the same here and use effective theory methods [ ]23[ ] to systematically expand in both 
small parameters, the coupling and {p^ — mf)/mj. 

The method has been worked out in detail and its viability has been shown in the case of a toy 
model with a single resonant heavy scalar particle [|^]. A first step towards considering pair pro- 



duction near threshold, albeit for the case of W bosons, has been described in Ref. [25]. However, 
the main features are the same as for top quark pair production. A consistent description requires 
the introduction of additional operators and modes. The incoming electrons and the decay prod- 
ucts of the top quarks are to be described by coUinear modes familiar from soft-coUinear effective 
theory ||2^]. The production and decay of the top quarks is described by new operators (ee) (z^'V), 
where e and e are collinear electron and positron fields respectively, and Y^W^b and xW^b, where 
W and b describe the decay products of the top quark. Electroweak corrections are split into various 
contributions. First there are the hard corrections to be taken into account by computing the match- 
ing coefficients of these operators to the desired accuracy. This corresponds to the computation 
of the matching coefficients in NRQCD and reproduces the so-called factorizable corrections. We 
are then left with matrix element corrections due to the still dynamical soft, potential and collinear 
modes which can be associated with the non-factorizable corrections. 

While there is still a long way to go to actually perform an explicit calculation, we are confident 
that this will be the method of choice to include the finite width effects in a systematic way. 

5. Conclusions 

There has been a lot of progress in the theoretical evaluation of the tt threshold cross section 
and the situation concerning the position of the peak of the cross section, relevant for the mass 
measurement is very well under control. A mass measurement with an error of 100-150 MeV is 
certainly achievable, as long as a threshold mass definition is used. Regarding the normalization 
of the cross section, relevant for a measurement of the width of the top quark and eventually its 

8 
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Yukawa coupling, the resummation of the logarithms has also substantially improved the situation. 
However, this will need further improvements to fully exploit the potentially very precise experi- 
mental data. Apart from the completion of the NNLL evaluation, the effects due to the finite width 
will have to be fully understood. Obviously, any progress towards a NNNLO evaluation is most 
welcome. Finally, let me conclude with an obvious but often forgotten remark: exploiting this 
unique opportunity to measure some Standard Model parameters to an unprecedented accuracy 
actually does require a dedicated experimental effort in that a future International Linear Collider 
will have to be run at the top threshold for a reasonable amount of time. 
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